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THE RESPONSE OF THE UPPER ATMOSPHERE TO 


PERTURBATIONS FROM DIFFUSIVE EQUILIBRIUM 
P. W. Blum 
ABSTRACT 

It is generally assumed that in the atmosphere above 120 km the deviations 
from diffusive equilibrium are small, though a minor constituent may show 
some deviation. The response of such a constituent to perturbations from dif- 
fusive equilibrium is qualitatively analyzed. It is shown that the magnitude of 
these deviations is mainly determined by the characteristic diffusion time. For 
a time-dependent perturbation the ratio of the characteristic time of the per- 
turbation to the characteristics diffusion time is an important parameter. 
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THE RESPONSE OF THE UPPER ATMOSPHERE TO 
PERTURBATIONS FROM DIFFUSIVE EQUILIBRRJM 

INTRODUCTION 

The atmosphere above a height of 120 km is generally considered to be 
close to diffusive equilibrium. Most atmospheric models are based on this 
assumption. Especially, the very good description of atmospheric densities 
given by the semi-empiric Jacchia model (Jacchia, 1971) that is based on dif- 
fusive equilibrium, shows that diffusive equilibrium of the atmospheric consti- 
tuents, save possibly hydrogen and helium, is a close approximation to the true 
state of the atmosphere. On the other hand, some of the more detailed atmos- 
pheric structure may be related to deviations from diffusive equilibrium (Mayr 
and Volland, 1971). 

The following processes must be considered as possible sources of devia- 
tions from diffusive equilibrium; 

1. Eddy diffusion 

2. Dissociation and recombination of oxygen 

3. Horizontal flow having a non-zero divergence. 

In this paper we shall obtain qualitative results concerning the response 
of a minor constituent to a perturbation from diffusive equilibrium. While for 
a purely qualitative analysis the nature of the soxirce of the deviation has no 
particular importance, it should be remarked that for the atmosphere above 
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120 km both eddy diffusion and recombination of atomic oxygen may be disre- 
garded with little loss of accuracy. 

For the recombination of atomic oxygen this is immediately apparent by a 
consideration of the cross section of the process. Regarding eddy diffusion we 
have to compare the eddy diffusion coefficient with the molecular diffusion co- 
efficient. The ratio of the eddy diffusion coefficient to the molecular diffusion 
coefficient is at 120 km 0.1 and at 160 km 0.01, if a height independent eddy dif- 
fusion coefficient is assumed . However, this height- independency cannot be 
maintained (ShimazaM, 1971), and a rather sharp decrease by one order of 
magnitude near a height of 105 km is probable. Taking this behaviour into ac- 
count, no eddy diffusion above 120 km, or even a slightly lower altitude, need 
be considered. 

In this paper we are only concerned with qualitative results, therefore we 
may use a simplified atmospheric model. Its properties are; 

1. The major constituent is in diffusive equilibrium and at rest. 

2. The minor constituent is in diffusive equilibrium at the upper boundary. 

3. The temperature profile is given. 


BASIC EQUATIONS AND DEFINITIONS 


Chapman and Cowling's diffusion equation for the vertical diffusion velocity 
v^ of two atmospheric constituents reads in their original notation 


''d = - Vj = - 
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Here v. are the constituent velocities, either absolute or relative to the center 
of mass, Di 2 the binary diffusion coefficient, n j the number densities and n their 
sum, n^^ the ratio n^/n, m. the molecular masses, p the density, p the total 
pressure and F. the external forces. For the constituents of the neutral 
atmosphere the latter are only the gravitational attraction and therefore equal. 

The first approximation of the diffusion coefficient given by Chap- 

man and Cowling is 
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where T is the temperature and a the effective collision diameter. A rear- 
rangement of (1) yields the form 


-L ^ _L _L _ /!JL ^ J_ _L 

rij 3z Dt \U 2 3z Dt // 


(3) 


where v^ is the center of mass motion and the scale heights. When the in- 
ertial forces are neglected compared to g and use is made of the assumption 
that the major constituent is in diffusive equilibrium, the simpler form 
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(4) 


results. 

For a constituent that is in diffusive equilibrium the expression in brackets 
of (4) vanishes and accordingly the diffusion velocities are zero. 
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The value of the bracketed expression for a constituent that is not in dif- 
fusive equilibrium will indicate the degree of its deviation from diffusive equili- 
brium. A dimensionless quantity d that will be a measure of the deviation will 
be defined by 


d 



3n 

Bz 
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^ H ^ T Bzj 


( 5 ) 


where the indices were dropped for simpler notation. 


CONSIDERATIONS OF DIMENSIONS 

Already from considerations of dimensions we may deduce possible ex- 
pressions for the characteristic time scales in the atmosphere. The quantities 
that appear in atmospheric dynamics are for a single constituent atmosphere 
k, T, H, m, g, n and the collision diameter cr^. From these quantities only few 
expressions having the dimension of time may be formed. If k, T or m appear 
in such a combination then they must appear as kT/m because in the result 
neither temperature nor mass may appear. But kT/m equals Hg and therefore 
we may limit ourselves to combinations of H, g, n and cr^. Of these only g in- 
cludes a time component, it must therefore be in the combination with an ex- 
ponent of -1/2. As exponents higher than 3 or 4 for the other quantities are 
improbable, only a rather narrow choice of combinations that have the dimen- 
sion of time is available. Two such combinations suggest themselves: ^ H/g 
and (/llg n)"^ . The first will be the characteristic time required for a gas 
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to return to barometric equilibrium if a large scale deviation from equilibrium 
exists, the second wili be the characteristic life-time of microscopic deviations 
that are equalized by molecular collisions. 

In diffusion processes two atmpspheric constituents take part and we have 
to add D ^2 to toe other quantities in order to obtain a time-like combination. 

The ratio of the molecular masses, i.e. the scale heights, may also appear as a 
dimensionless coefficient. the dimension of cm^ sec“^ and is in itself 

a combination of H, g, n and cr^. Again no unique determination of the character- 
istic time from considerations of dimensions seems possible, but H Vd suggests 
itself as the simplest physically plausible form. Here the gravitational force, 
the temperature and the collision diameter have all been absorbed in • 

EQUALIZATION OF AN INITIAL PERTURBATION 

FROM DIFFUSIVE EQUILIBRIUM 

In the following we shall use the subscript m for the major constituent, i 
for the minor constituent, o for a time independent component, b for the initial 
distribution and primes to indicate partial differentiation with respect to the 
vertical coordinate. 

For the major constituent we have 


n 


n 


III m 3T n 

n + — + — = 0 

H T 3z 


( 6 ) 


For the minor constituent we assume an initial deviation from diffusive equili- 
brium. In the absence of a continuous perturbation the minor constituent 
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will finally return to diffusive equilibrium and have then the density distribution 
n^ where 


n 



d ^ 

T 3 


( 7 ) 


The density of the minor constituent may be written as 


n(z, t) = rij (z) + (z) 


(8) 


where n^ is the part of the initial distribution that is not in diffusive equilibrium. 
Ub may have positive or negative sign. As we are only concerned with qualita- 
tive resvdts we shall assume the atmosphere to be isothermal. This simplifi- 
cation will not affect the order of magnitude of the results, but will simplify 
the analytic treatment. 

Perturbations that have a vertical extension of one scale height will be 
considered. We have therefore -z^ % H where z^ is the upper boimdary and 
Zq the lower boundary. At the upper boundary z^^ the minor constituents will 
always be in diffusive equilibrium and have a density equal to the final density. In 
order to realize this model we define an effective scale height 

H^=H + AH = H ^1 + = H (1 + / 3 ) (9) 

with the properties 
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ah (z^) = 0 

fi{z j) = 0 


( 10 ) 



n'+-?-=0 (12) 

H 

e 

This definition assures the gradual disappearance of the deviation with time. 

We obtain the partial time-derivative and the diffusion velocity 



Substitution in the continuity equation 


( 13 ) 
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gives 
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, n in nfi 
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therefore 


I H / H H H 


(18) 


The last term of (18) may be neglected as it is the product of small quantities. 
From (16) 


-a n. = — D 

3z \ H H y 


(19) 


We differentiate each term of the product on the right hand side of (19) and sub- 
stitute after the differentiation the boundary values, i.e. z-z^ = -H. From (2) it 
is evident that 

3D _ D 

3z H (20) 
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This yields 


an. = - — — n ,Pn + + 

b jjji d^O ^2 u2 


(21) 


H' 


All three terms are of the order of magnitude D n ^ The sum of the first 

two terms may introduce a numerical factor of the order of unity in our result. 
Disregarding this factor we have 
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( 22 ) 
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b 


The positive sign of a is assured by the negative sign of/Sg/n^,. The magnitude 
of the coefficient nj/^o/n^ is determined from the solution of (7) and (12) where 
the integration constants are adjusted to yield equal densities at the upper 
boundary. Using this relation the coefficient becomes 


_ ^0 
rij, 1 - exp 


(23) 


and is therefore close to unity for small 

We may according to the above estimate define as the characteristic time 
of diffusion in the upper atmosphere t ^ = H^/D sec. 

We have obtained the result that after H^/D sec the deviation from diffusive 
equilibrium has fallen to 1/e of its initial value. 

Some confusion may arise from the proportionality of the diffusion velocities 
to the diffusion coefficient as expressed in the diffusion equation (1). This seems 
to indicate an increase of diffusion velocities with height. This is only true for 
the initial diffusion velocities that arise from a instantaneous perturbation. As 
the diffusion velocities decrease exponentially with time with an exponent -D/H ^ 
after a time span of the order of the diffusion time, the velocities at the higher 
altitude will be lower than at the lower altitude. 
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Figure 1 shows the height profile of D and t^ as derived from the Jacchia 
model for To, = 1000° and atomic oxygen and nitrogen as the constituents. 

CONTINUOUS PERTURBATIONS FROM DIFFUSIVE EQUILIBRIUM 

Deviations from diffusive equilibrium that arise from a continuous source 
will now be considered. Such sources of deviation are the horizontal divergence 
of the horizontal particle flux or the production and loss term of particles. First 
the case of time-independent perturbations will be treated. This corresponds 
to the diurnally averted deviations caused by the diurnal average of the above 
sources. While the diurnal average of the horizontal divergence is probably 
small compared to its instantaneous values, this is not the case for the produc- 
tion of atomic oxygen which has an average value comparable to its instantaneous 
values as recombination above 120 km is negligible. According to Shimazaki's 
results the production and loss term is of the order of 10~® n sec~^ cm~^ in the 
thermosphere above 120 km. We shall always assume that the perturbation is 
proportional to the number density, an assumption that is self-evident for the 
horizontal divergence as well as the production and loss term. 

We shall assume the perturbation to be F(z) n sec"^ cm"^, where F(z) will 
be only slightly height dependent. 

The continuity equation for the minor constituent for the time- independent 
mode becomes 

^ (vjn)o = F(z) no (24) 
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or 
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--f 


F(z) no dz = FH 


(25) 


where F is a height average of F(z) over one scale height. The deviation d is 
given by (5). Accordingly 


d = 


VdH 

D 


FH^ 

D 


(26) 


We have therefore the result that the magnitude of the deviation from dif- 
fusive equilibrium is proportional to the characteristic diffusion time and the 
magnitude of the perturbation. 

TIME DEPENDENT PERTURBATIONS 

In the case of time dependent perturbations the deviation from diffusive 
equilibrium will be dependent on the characteristic diffusion time and the 
characteristic time of the perturbation. We shall assume a perturbation of the 
form 


F^(z, t) = F(z) (27) 

The density will be given by 

n = Oq (1 + d (28) 
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A similar treatment as in the time-independent case yields 


(v<,n)=F(z)ne‘"' 


(29) 
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because 


■^0 “ 
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Integration of (30) resvdts in 


lO) 


r*i , _ 

I doQ dz - D n^d = F ri(jH 


(33) 


where the right hand side was linearized as d is assumed to be small. 


icoH _ _ FH 
D D 


(34) 


We define a ratio r of the characteristic times as 



p 


and obtain 


(35) 
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(36) 


d' - 


i rd 
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If r « 1 the equation reduces to 


d' =-Il± (37) 

H 

and integration over a scale height yields 

d = Ft^ (38) 

This expresses again that the deviation is proportional to the amplitude of the 
perturbation and the characteristic diffusion time. 

In case the ratio r has a magnitude that approaches, or exceeds, unity, we 
may not neglect the second term of (36). In this case the solution is given by 


exp - 


(- fi-) r 


F(z)t^ 

H 


exp 



(39) 


where we have made use of the vanishing deviation at the upper boundary. In this 
case, which is applicable to the height region below about 160 km, the resulting 
deviation depends on the ratio of the characteristic times and- the characteristic 
diffusion time. 

We may solve (39) formally by observing that both r and t^ decrease ex- 
ponentially with altitude. The result will be a deviation that has an altitude 
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dependent phase factor and an amplitude that is proportional to the perturbation 
amplitude F(z) and the diffusion time H VD. 


d » Ft^ 


Ft^ 


(40) 


RESULTS 

1. The upper atmosphere responds to perturbations from diffusive equili- 
brium with a character time H /D. 

2. Perturbations that extends over a height range of the order of a scale 
height are damped with a factor exp (-D/H ^ • t), 

3. Deviations from diffusive equilibrium that arise due to continuously act- 
ing sources are proportional to the amplitude of the perturbation and the charac- 
teristic diffusion time. 

4. In the height region where the characteristic time of the perturbation is 
of the order of the characteristic diffusion time, or of smaller order, a height 
dependent phase shift between the deviation and the source is introduced. 
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